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ARTICLE INFO ABSTRACT 



The 5386 nucleotide bacteriophage 0X1 74 genome has a complicated architecture that encodes 11 gene 
products via overlapping protein coding sequences spanning multiple reading frames. We designed 
a 6302 nucleotide synthetic surrogate, 0X174.1, that fully separates all primary phage protein coding 
sequences along with cognate translation control elements. To specify 0X174.1f, a decompressed 
genome the same length as wild type, we truncated the gene F coding sequence. We synthesized DNA 
encoding fragments of 0X174.1f and used a combination of in vitro- and yeast-based assembly to 
produce yeast vectors encoding natural or designer bacteriophage genomes. We isolated clonal 
preparations of yeast plasmid DNA and transfected E. coli C strains. We recovered viable 0X174 
particles containing the 0X174.1f genome from E. coli C strains that independently express full-length 
gene F. We expect that yeast can serve as a genomic 'drydock' within which to maintain and 
manipulate clonal lineages of other obligate lytic phage. 

© 2012 Elsevier Inc. All rights reserved. 
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Introduction 

Small icosahedral virus particles containing single stranded 
DNA genomes encoding various Microviridae have been found in 
many environments, from seawater to the human gut (Roux et al., 
2012). Given its intrinsic natural importance plus accumulated 
relevance as an experimentally tractable model system, bacter- 
iophage 0XI 74 is now the best-studied member of the Microvir- 
idae (Hayashi et al., 1988; Fane et al, 2005). 0X174 is also widely 
used to explore topics ranging from adaptive evolution, to 
molecular self-assembly, molecular systems biology, and virus 
particle transmission (Wichman and Brown, 2010; Gordon et al., 
2012; Markely and Yin, 2007; Julian et al., 2010). 

One notable feature of the 0XI 74 genome is the extent to which 
the coding sequences for seemingly independent gene products 
overlap (Sanger et al., 1977). For example, gene E, whose encoded 
product is responsible for host cell lysis, is encoded entirely within an 
alternate reading frame internal to gene D, which encodes an external 
scaffold protein (Barrell et al., 1976; Fane et al., 2005). Overall, 16.8% 
of the 0XI 74 genome encodes for more than one gene, and there are 
two genomic regions in which three genes overlap (Fig. 1A, top). 

For broader context, ~30% of prokaryotic genes overlap with the 
coding sequences of neighboring genes (Kingsford et al., 2007), 
and several models have been proposed to explain the existence of 
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overlapping genetic elements. For example, Yokoo and Oshima (1979) 
hypothesized that the 0X174 genome encodes information from an 
advanced alien civilization, with overlapping genes as the most likely 
candidate regions in which to find such messages. Researchers have 
also considered co-regulation of gene expression, genetic coding effi- 
ciency as driven by genome replication or packaging constraints, and 
coupling of evolutionary selection(s) as terrestrial factors leading to 
overlapping elements (Chirico et al., 2010). However, direct experi- 
mental tests demonstrating the spontaneous emergence of overlap- 
ping genes and other genetic elements have been difficult to imple- 
ment given the starting complexity of already evolved natural systems. 

Previously, using the bacteriophage T7 genome as a test case 
we introduced and explored the idea that natural genetic systems 
might be 'refactored' in order to produce engineered surrogates 
that are easier to study and apply (Chan et al., 2005). Our 
refactoring process reflected six engineering goals: (1) well- 
defined boundaries for all genetic elements, (2) no overlapping 
elements, (3) only one encoded function per element, (4) enable 
precise and independent manipulation of elements, (5) a genome 
that is possible to construct, (6) a genome that encodes a viable 
system. We constructed 12,179 base pairs of the T7.1 genome and 
recovered viable chimeric phage, albeit with significantly reduced 
fitness. Subsequent phage evolution experiments by others 
reported that ~35% of the T7.1 design changes were lost while 
selecting for mutant phage that recovered to a near wild-type 
fitness level (Springman et al., 2012). Taken together, these two 
studies suggest that genome refactoring to remove overlapping 
elements followed by experiments that attempt to select for 
"repacking" might be fruitful. 
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Fig. 1. Design of a fully decompressed bacteriophage 0X174 genome. (A) Reading frame comparison of the wild type and decompressed oX174.1f. Wild-type genome reading 
frames +1 to +3 are given from top to bottom. (B) Expanded view of the wild-type 0X174 genes A, K, C junction. (C) Decompression of genes A, K, and C. To separate gene K from 
gene A (top) the start codon of K was changed from ATG to ACG and its RBS weakened. To decompress gene C from gene A (top) the start codon was changed from ATG to GTG and 
its RBS weakened. To decompress gene K from gene B (middle) the start codon was changed from ATG to ATT and its RBS weakened. To decompress gene C from gene K (bottom) 
the start codon was changed from ATG to ACG and its RBS weakened. Start and stop codons are green and red, respectively, while overlapping start/stop codons are blue. RBSs are 
underlined. Changes made during gene decompression are shown as capital letters. 



However, given our prior experiences with T7, we first sought to 
develop a new system and methods that would be better suited for 
developing and applying ideas related to refactoring to obligate lytic 
phage. For example, many DNA fragments encoding phage genetic 
elements cannot be readily propagated using typical E. coli cloning 
strains and vectors. As a second example, we were limited in working 
on T7.1 by concerns that any individual genetic design change might 
produce a non-viable genome that could not be physically recovered, 
propagated, or amplified; since refactoring projects often attempt to 
change hundreds of nucleotides simultaneously (Temme et al, 2012), 
non-viable designs that cannot be independently cloned or propa- 
gated are difficult if not impossible to work with. 

There have been notable advances in the synthesis, assembly, and 
handling of nucleic acids over the past decade (Carr and Church, 
2009). For example, the genomes of several viruses and one microbe 
have been assembled from chemically synthesized oligonucleotides, 
with viable systems being recovered in each case (Cello et al, 2002; 
Smith et al, 2003; Gibson et al, 2010; Yang et al, 2011). Also of note, 
a bacterial genome was assembled via transformation-associated 
recombination in yeast, wherein the yeast-based bacterial genome 
might be maintained and further manipulated without immediate 
concern for the viability of the microbe resulting from the 
so-encoded bacterial genome (Gibson et al, 2008). 

Thus, we sought to establish methods enabling the synthesis, 
assembly, and recovery of a bacteriophage genome via yeast. We 
first choose to work with bacteriophage 0X174 given its small 
circular genome. We were then further motivated to explore if all 
genetic element overlaps within the 0X174 genome might be 
eliminated and shown to be non-essential, generating the first 
fully decompressed phage genome. 



Results and discussion 

Designing a decompressed bacteriophage 0X2 74 

We designed a completely decompressed version of 0X174 in 
order to test the limits of refactoring with a highly overlapping 
natural genetic system. Our first design, 0X174.1, eliminated all gene 
overlaps while preserving the use of native ribosome binding sites 
(RBSs) and promoters as much as possible (Fig. 1A, bottom). Our 
genetic decompression algorithm was similar to that used with T7: 
(i) overlapping coding sequences were fully separated, with sequence 
duplication as needed, into physically distinct open reading frames 
and (ii) vestiges of regulatory elements such as start codons and 
legacy RBSs were weakened or erased via silent point mutations so 
as to make translation initiation of legacy reading frames less likely 
(Fig. IB and C). We generated synonymous mutations to remove 
remnant RBSs by reducing the frequency of A and G nucleotides. We 
used a published RBS strength calculator to confirm that the 
predicted strengths of so-changed RBSs decreased from wild-type 
strengths (Salis et al, 2009). Again following Chan et al. (2005), in 
making synonymous changes we sought to maintain codon usage 
bias but, when required, choose higher frequency codons. 

Next, we used these design rules to separate genes B, K, C from 
each other and gene A, gene D from gene C, and genes E and J from 
gene D (Fig. 1, Supplementary Files). We also implemented 22 
simultaneous point mutations from the wild-type sequence in order 
to disrupt legacy start codons and RBSs (Table 1 ). We referred to the 
so-changed coding and RBS sequences via a "letter.l" notation; for 
example, gene Al replaces gene A. In total, the process of genome 
decompression added 909 nucleotides to the wild-type genome. 



280 P.R. Jaschke et al. / Virology 434 (2012) 278-284 



Table 1 

Single nucleotide changes in the decompressed oX174.1f design to remove legacy gene sequences. 



Position a 


Refactored 


Refactored 


Position 11 


Original 


Original 


Effect 




gene 


sequence 




gene c 


sequence 




1 149 


A.l 


T 


5066 


B 


A 


Gene B RBS remnant weakened 


1152 


A.l 


c 


5069 


B 


T 


Gene B RBS remnant weakened 


1158 


A.l 


c 


5075 


B 


A 


Gene B remnant start codon changed ATG- > GTG 


1509 


A.l 


T 


40 


K 


G 


Gene K RBS remnant weakened 


1512 


A.l 


c 


43 


K 


T 


Gene K RBS remnant weakened 


1521 


A.l 


c 


52 


K 


T 


Gene K remnant start codon changed ATG- > ACG 


1590 


A.l 


c 


121 


c 


T 


Gene C RBS remnant weakened 


1593 


A.l 


c 


124 


c 


T 


Gene C RBS remnant weakened 


1602 


A.l 


c 


133 


c 


A 


Gene C remnant start codon changed ATG- > GTG 


1974 


B.l 


T 


42 


K 


A 


Gene K RBS remnant weakened 


1985 


K.1 RBS 


T 


53 


K 


G 


Gene K remnant start codon changed ATG- > ATT 


2069 


K.1 


c 


122 


c 


G 


Gene C RBS remnant weakened 


2072 


K.l 


C 


125 


c 


G 


Gene C RBS remnant weakened 


2081 


K.1 


C 


134 


c 


T 


Gene C remnant start codon changed ATG- > ACG 


2428 


C.l 


C 


378 


D 


A 


Gene D RBS remnant weakened 


2431 


C.l 


c 


381 


D 


T 


Gene D RBS remnant weakened 


2440 


C.l 


T 


390 


D 


A 


Gene D remnant start codon changed ATG- > TTG 


2626 


D.l 


A 


557 


E 


G 


Gene E RBS remnant weakened 


2638 


D.l 


C 


569 


E 


T 


Gene E remnant start codon changed ATG- > ACG 


2905 


D.l 


C 


836 


J 


A 


Gene J RBS remnant weakened 


2908 


D.l 


T 


839 


J 


A 


Gene J RBS remnant weakened 


2919 


E.l RBS 


T 


850 


J 


G 


Gene J remnant start codon changed ATG- > ATT 


3378-3385 


After J/F-intergenic 


GCGATCGC 








AsiSI site added J/F-intergenic 



a Coordinates of refactored oX174.1f (see Supplemental Data for sequence file), beginning at the start of the H/A-intergenic region. 
b Coordinates according to original 0X1 74 sequence numbering scheme. 

c Sequence belonging to original gene RBS or start codon although it is embedded in another gene in the wild type sequence. 



We next addressed practical constraints arising from the length of 
DNA that can be physically packaged within a 0XI 74 capsid without 
impacts to reproductive fitness. Previous work has shown that the 
length of a 0X174 genome, when packaged in vitro, must be kept 
within a few percent of the 5386 nucleotide wild-type length in order 
to avoid any significant fitness decrease (Aoyama and Hayashi, 1985). 
Similar results were shown in vivo (Russell and Muller, 1984). To 
reduce the decompressed genome length we removed the first 916 
nucleotides of gene F, encoding the coat protein (Air et al„ 1978). We 
chose gene F because a plasmid containing a restriction fragment 
encoding wild-type gene F was able to complement two conditional 
gene F mutations (Avoort et al„ 1983). Additionally, the gene F coding 
sequence is greater than the total of the combined increases needed 
to implement the 0X174.1 genome design. The truncated gene F 
version of the decompressed genome was named oX174.1f. 
To complement 0X174. If when transformed into host cells we 
designed a medium copy vector expressing gene F under control 
of a rhamnose-inducible promoter (Fig. SI). 

Building bacteriophage genomes in yeast 

Our workflow to construct bacteriophage genomes in yeast 
consisted of several steps (Fig. 2). After designing a genome on a 
computer we obtained all genomic elements to be assembled via 
DNA synthesis or PCR cloning. We then used PCR amplification 
to add homologous overhangs to adjacent elements and any short 
functional sequences, such as restriction enzyme recognition 
sites, required for later steps. We mixed all DNA fragments 
together with a yeast/£. coli shuttle vector backbone and, if 
needed, performed in vitro isothermal assembly (Gibson et al., 
2009) followed by transformation into yeast spheroplasts 
(Kouprina and Larionov, 2008; Orr-Weaver et al., 1983). We 
screened yeast colonies using multiplex PCR followed by plasmid 
purification from positive clones. These plasmids were subse- 
quently digested to release the phage genome and ligated to re- 
form a circular phage genome. Re-circularized genomes were 



then transfected into suitable host cells. Finally, when our design 
and assembly process resulted in viable phage, we picked plaques 
and sequenced genomic DNA from clonal lineages. 

Transplantation of wild-type 0X174 to and from yeast 

To test the validity of our assembly workflow and to determine 
whether bacteriophage genomes propagated in yeast remain 
viable, we first attempted to capture and then recover a wild- 
type 0X174 genome from yeast. We used PCR to linearize the 
yeast shuttle vector pRS425 while also adding overhangs homo- 
logous to the ends of a linearized wild-type 0X174 genome plus 
AsiSI restriction sites (Fig. 2A). Both DNA fragments were co- 
transformed into yeast spheroplasts. Plasmids from positive 
screened clones (Fig. 2B) were processed via restriction digestion 
and ligation and transformed into E. coli C to establish that a 
recovered 0XI 74 genome would still initiate infection. We recov- 
ered equivalent size plaques following overnight plating for both 
a wild-type positive control and also with 0X174( AsiSI) (Fig. 2C, 
Table 2). We did not recover plaques from transformations for 
which the yeast vector backbone was not removed via AsiSI 
digestion. Overall, of the 16 yeast clones we screened, 11 were 
positive by PCR for both the phage-vector junctions and 7 of these 
clones produced phage, giving an overall yeast assembly followed 
by redeployment process efficiency of ~44%. 

Assembly and recovery of a decompressed bacteriophage oXl 74 

We next sought to construct and test a 0X174.1f genome. We 
assembled the oX174.1f genome from nine DNA fragments encod- 
ing independent phage elements plus the shuttle vector backbone. 
Because assembly efficiencies decline with increasing numbers of 
fragments (Liang et al., 2012), we used longer 30 base pair 
homologous fragment overhangs. Yeast-based assembly alone 
yielded no functional phage. Therefore, to increase the proportion 
of complete assemblies, we used an in vitro isothermal assembly 
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Fig. 2. Building bacteriophage with yeast. Complete workflow for assembling and recovering a bacteriophage genome from a combination of natural and synthesized DNA 
parts. From left to right, the DNA sequence of a new bacteriophage genome is designed using software tools. The required parts are sourced from either natural sequences or 
DNA synthesis. Parts are amplified by PCR. Assemblies involving more than two parts are assembled in vitro followed by transformation into yeast spheroplasts. An auxotrophic 
marker encoded on the yeast vector is selected for. Transformants are screened using colony multiplex PCR with primer sets spanning key assembly junctions. Cultures of 
positively screened colonies are grown and plasmid DNA is purified directly from cell pellets. Phage genomic DNA is then released from the yeast vector, religated, and 
transformed into bacterial host cells. Resultant plaques are sequenced to verify correctness of the assembly. (A) A linearized wild-type oXl 74 genome assembled with yeast 
vector pRS425. Flanking AsiSI restriction sites were added to the assembly via PCR primers used to linearize the vector, so that the phage genome can be released from the 
vector and re-circularized prior to bacterial transfection. (B) Multiplex PCR screening of pRS425(oX174( AsiSI)) assemblies. The left schematic shows amplicon locations 
for junction testing. Letters (right) note expected amplicon sizes for correct assemblies. (C) Infection of E. coli C with 0X174(AsiSI) genomes recovered from yeast. 



Table 2 

Plaque sizes of wild type and refactored oX174.1f strains grown with varying 
amounts of gene F induction. 



Host strain Phage Rhamnose (uM) Plaque size' 1 



E. coli C 


Wt 0X174 




6(1) 


E. coli C 


0l74X( AsiSI) 




6(1) 


E. coli C 


Refactored 0X174.1f 




n/a 


E. coli C(p804(Gene F)) 


Wt 0X174 


0 


4.9(0.7) 


£ .coli C(p804(Gene F)) 


Wt 0X174 


1 


4.6(0.6) 


E. coli C(p804(Gene F)) 


Wt 0X174 


10 


4(1) 


E. coli C(p804(Gene F)) 


Wt 0X174 


100 


4(1) 


E. coli C(p804(Gene F)) 


Wt 0X174 


1000 


4.4(0.8) 


E. coli C(p804(Gene F)) 


Refactored 0X174.1f 


0 


n/a 


£. coli C(p804(Gene F)) 


Refactored 0X174.1f 


1 


0.6(0.1) 


E. coli C(p804(Gene F)) 


Refactored 0X174.1f 


10 


0.9(0.2) 


E. coli C(p804(Gene F)) 


Refactored oX174.1f 


100 


2.2(0.7) 


£ .coli C(p804(Gene F)) 


Refactored oX174.1f 


1000 


4.1(0.6) 



a Plaque sizes in mm + standard deviation in brackets. 



step (Gibson et al„ 2009) prior to yeast transformation. We 
screened resulting yeast clones using two multiplex primer pairs 
targeting both phage-vector junctions (Fig. 3A). Positive clones 
were processed as before and transformed into an E. coli C strain 
expressing gene F under the control of rhamnose induction. Full 
genome sequencing of phage DNA recovered from a clonal lineage 
cultured from a plaque confirmed that all designed changes to the 
wild-type genome were present in 0X174. If. Additionally, in gene 
G, where our design recapitulated the wild-type sequence (Sanger 
et al., 1977) we identified two unexpected nucleotide changes: a 
GAT (Asn) to AAT (Asp) substitution in the 1 13th codon, and an AAC 
to AAT silent mutation in the 139th codon. Both of these differences 



from our design likely arose from sourcing of the DNA fragment 
encoding gene G via PCR from our laboratory phage stock (as 
opposed to de novo DNA synthesis). Both polymorphisms are 
consistent with observed variation in the 0X174 genome, as 
Smith et al. (2003) also reported these exact two changes from 
sequencing of laboratory stocks of 0X174. 

Initial phenotypic characterization of 0XI 74. If 

We checked for first-order phenotypic effects, if any, resulting 
from the 0X174.1f genome architecture. For reference, separation 
of bacteriophage T7 early genes resulted in a 50% reduction in 
plaque size and lysis times (Chan et al., 2005). First, we compared 
wild-type and decompressed phage plaque formation on an E. coli 
C host expressing various amounts of gene F under the control of 
a rhamnose inducible promoter (Fig. 3B). Wild-type plaque sizes 
did not change significantly with varying gene F induction 
(Table 2). No plaques were observed for 0X174.1f when 0 urn 
rhamnose was added. 0X174.1f plaque sizes increased from ~ 7.5- 
fold smaller than wild-type at 1 urn rhamnose to virtually 
identical to wild-type at 1000 p.m (Table 2). These results suggest 
that gene F expression levels are closely coupled to plaque size. 

We also assayed the growth properties of oX174.1f in liquid 
culture. As expected, the decompressed 0X174.1f phage was unable 
to lyse E. coli C lacking the gene F plasmid (Fig. 3C, left panel). We 
observed that both wild-type 0X174 and 0X174.1f phage lysed 
E. coli C (pJ804(Gene F)) with similar efficiency (Fig. 3C, right panel). 

Taken together, when sufficient gene F is supplied in trans the 
first-order phenotype of 0X174.1f is indistinguishable from the wild 
type. In an artificial fashion, our decompression of the 0XI 74 genome 
has allowed us to travel 'backwards in time' regarding Microviridae 
evolution. For example, previous bioinformatic analysis has suggested 
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Fig. 3. Assembly and characterization of a the fully decompressed oX1 74.1 f genome. (A) Multiplex PCR screening of pRS425(oXl 74.1 f) assemblies. Schematic on left shows 
locations of PCR amplicons. Letters on right denote expected amplicon sizes for correct assembly of the first and last phage parts with the yeast vector. (B) Characterization 
of wild-type 0X1 74 and decompressed 0XI 74.1 f plaque size in £. coli C strain expressing gene F under the control of a rhamnose inducible promoter. (C) Characterization of 
wild-type 0X174 and decompressed 0X1 74.1 f liquid lysis dynamics. Left panel depicts infection of E. coli C; right panel depict infection of E. coli C(pJ804(Gene F)) in media 
supplemented with 10 uM rhamnose. Cells were grown in TK broth with 5 mM CaCl 2 and 10 mM MgCl 2 added just prior to infection. Assays were carried out at 37C with 
an initial average MOI of ~5. Uninfected control (red), decompressed oX174.1f (light blue), and wild-type 0X174 phage (dark blue). Error bars give standard deviations 
from four replicates. 



that 0X174 genes evolved in a step-wise fashion, whereby gene A 
first expanded its 3' end to overlap gene B, followed by overprinting 
of genes K and E (Pavesi, 2006). Such findings are supported by the 
non-essential nature of genes A*, K and E (Colasanti and Denhardt, 
1987; Gillam et al., 1985; Pollock et al, 1978). Our aggregate finding 
that no essential information is lost by full oX174 gene decompres- 
sion may further support the hypothesis that genome compression 
results from evolutionary pressure to package more genetic informa- 
tion within a limited capsid size (Chirico et al., 2010). Now that the 
0X174 genetic sequences have been separated from one another 
further tests of such evolutionary relationships can be pursued via 
direct experiments. 



Materials and methods 

Bacterial, yeast and bacteriophage strains and propagation 

Wild-type 0XI 74 (ATCC # 13706-B1 ) and E. coli C host (ATCC # 
13706) were used throughout. All competent cells were made 
by the TSS method (Chung et al., 1989). Bacterial cells were 
propagated with LB and TK broth (Fane and Hayashi, 1991 ). E. coli 
C (pJ804(Gene F)) was routinely grown with carbenicillin at 
50 u.g/mL concentration. Yeast strain W303 was used. 0XI 74 
was propagated following published methods (Fane and Hayashi, 
1991). 



DNA design and manipulation and sequencing 

0X174 wild-type sequence NC_001422.1 was used throughout 
as the reference 0X174 sequence. All wild-type sequences are 
numbered as in NC_001422.1 (from the single PstI cleavage site) 
as originally designated (Sanger et al., 1977). All refactored 
sequence coordinates are referenced from the first nucleotide of 
the H/A-intergenic region. A plasmid Editor (ApE, <http://biolo 
gylabs.utah.edu/jorgensen/wayned/ape/>) was used to visualize 
and manipulate all sequences. The Primer3Plus web interface 
( <http://primer3plus.com/cgi-bin/dev/primer3plus.cgi> ) (Rozen 
and Skaletsky, 2000), NCBI Primer-BLAST (Ye et al., 2012), and 
NetPrimer (Premier Biosoft) were used for designing and analyz- 
ing oligonucleotides. Synthetic DNA for genes A.1, B.l, K.l, C.l, 
D.l, E.l and the pJ804(Gene F) plasmid were synthesized by 
DNA2.0, Inc. The decompressed oX174.1f genome was sequenced 
via Sanger methods by Elim Biopharmaceuticals Inc., and we 
assembled sequencing reads into a contig using CodonCode 
Aligner. Average sequencing depth was greater than 3x at all 
positions. 

Building bacteriophage genomes in yeast 

To capture the 0X174 genome inside yeast we first linearized 
the yeast/E. coli shuttle vector pRS425 (Christianson et al., 1992) 
by using PCR amplification with opposite facing primers at 
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coordinates 2957 and 2971 Genbank Accession #1103452.1, leaving 
out a small part of the plasmid between the primers. The forward 
primer 1 (Table SI) encodes (5' -3') 20 nt of 0X174 homologous 
sequence, AsiSI site, and pRS425 homologous sequence. The reverse 
primer 2 was the same layout as the forward primer, both binding 
within lacZ. Dpnl digestion followed amplification. 

Since the native 0X174 genome is circular and single-stranded, 
we purified a dsDNA intermediate form of the genome (called RFI 
form) from infected cells (Burch et al., 1999). This purified circular 
dsDNA was linearized by PCR using opposite facing PCR primers 
(as for vector linearization above) that bind to the Gene H/A- 
intergenic region boundary. The forward primer 3 binds at the 
beginning of the H/A-intergenic region while the reverse primer 
4 binds at the 3'-end of Gene H. This site of linearization of the 
0X174 genome was chosen because the region of the genome: 
(1) does not encode any known genes, (2) is upstream of the P A - 
promoter and origin of viral strand replication located in Gene A 
(Hayashi et al., 1988), and (3) the junction is not known to encode 
any other function. 

Linear 0X174 genomes and pRS425 were co-transformed into 
yeast spheroplasts (strain W303) (Kouprina and Larionov, 2008; 
Orr- Weaver et al., 1983). Multiplex screening of the assemblies 
were performed on resuspended yeast colonies directly picked 
from 3 to 5 days re-grown spheroplast transformations. Yeast 
colonies were resuspended in 50 uL Zymolyase 20 T (MP Bio- 
chemicals) and 0.05 uL beta-mercaptoethanol, then incubated for 
1 h at 37 C. Primer sets for multiplex PCR were designed using 
muPlex (Rachlin et al., 2005). See Supplementary materials for 
primer sequences 5-8. Qiagen Multiplex PCR Plus kit was used for 
all screens. For 0X174( AsiSI) assembly screening the following 
PCR conditions were used: 95C/5:00, 35 cycles of 95C/0:30; 65C/ 
1:30; 72C/1:00, 68C/10:00, 8C/hold. 1-2 uL of yeast resuspen- 
sions were used for PCR template, while 1 ng of purified 0X1 74 
RFI (New England Biolabs) and pRS425 control templates used 
per 10 pi reaction. We cultured positive clones and purified 
plasmid DNA using a modified miniprep procedure (Singh and 
Weil, 2002). 

The decompressed oX174.1f genome was built identically to 
0XI 74( AsiSI), but with the following exceptions. Yeast vector 
pRS425 was linearized by Smal digestion, followed by gel pur- 
ification. The middle of gene F to the end of the H/A-intergenic 
region (coordinates 1917-3980) and gene J (coordinates 832- 
1000) parts of the assembly were amplified from wild-type 0X174 
RFI DNA generated from ATCC # 13706-B1. All other parts were 
amplified from plasmids synthesized by DNA2.0. All parts were 
assembled in vitro using the Gibson isothermal assembly method 
(Gibson et al., 2009), with 10 fmol of each part for 15 min at 50C, 
prior to yeast spheroplast transformation. Yeast clones were 
screened with multiplex primers 9-12. We initially screened 
yeast clones with denser multiplex PCR amplicons covering all 
part-part assembly junctions but dropped this expanded testing 
since we found the results to be difficult to interpret due to the 
highly repetitive nature of the 0X174.1f genome. Out of the 48 
yeast clones we screened, 18 were positive by PCR for both phage- 
vector junctions and only 1 of the prepared and processed vectors 
resulted in plaque formation, giving an overall success rate of 
~2%. We picked plaques, produced liquid culture phage lysates, 
and sequenced single-stranded genomic DNA from purified 
0XI 74.1 f phage particles. 

0X3 74 plaque assays 

For plaque assays, we added 100 ul phage in BE or 100 ul 
0X174 transformed E. coli C cells, plus 100 ul untransformed host 
E. coli C grown in LB broth, to 5 mL 0.7% TK agar at 46C. This 
mixture of cells and agar was poured over 1 5 mL 1 .2% TK agar 



plate and incubated at 37C for 12 h. For E. coli C(pJ804(Gene F)), 
antibiotics and rhamnose were added to the bottom agar in 
amounts needed for the total volume of bottom plus top agar. 

Lysis curves 

To obtain lysis curve data, 96-well plate experiments were set 
up in a Victor3 plate reader at 37C. A600 measurements were 
taken every 5 min, and MATLAB was used to process the data. The 
volume of each 200 ul well was comprised of phage, E. coli 
C(p804(Gene F)), 5 mM CaCl 2 , 10 mM MgCl 2 , variable amounts 
of rhamnose, and TK broth for the remaining volume. 
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